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SCIENCE ADVANCES | RESEARCH ARTICLE 


ECOLOGY 


Herbivores at the highest risk of extinction among 


mammals, birds, and reptiles 


Trisha B. Atwood'*", Shaley A. Valentine'*, Edd Hammill’, Douglas J. McCauley”, 
Elizabeth M. P. Madin’, Karen H. Beard, William D. Pearse®” 


As a result of their extensive home ranges and slow population growth rates, predators have often been per- 
ceived to suffer higher risks of extinction than other trophic groups. Our study challenges this extinction-risk 
paradigm by quantitatively comparing patterns of extinction risk across different trophic groups of mammals, 
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birds, and reptiles. We found that trophic level and body size were significant factors that influenced extinction 
risk in all taxa. At multiple spatial and temporal scales, herbivores, especially herbivorous reptiles and large-bodied 
herbivores, consistently have the highest proportions of threatened species. This observed elevated extinction 
risk for herbivores is ecologically consequential, given the important roles that herbivores are known to play in 


controlling ecosystem function. 


INTRODUCTION 

Over the past 500 years, at least 368 vertebrate species have gone ex- 
tinct (1), most as a result of direct (e.g., hunting) and indirect (e.g., 
habitat loss) human influences (2-4). Furthermore, 18% of extant 
vertebrates have been identified as threatened with extinction (1), 
and the expansion of anthropogenic disturbances combined with 
climate change is expected to exacerbate species declines (5). Be- 
cause species drive ecosystem functions, this species loss, like the 
historical mass extinction events before it, is likely to profoundly 
influence the ecology of our planet (6, 7). Historical extinctions have 
taught us that the disproportionate loss of large-bodied species and 
species within a specific trophic group (predators, herbivores, or 
omnivores) has substantial ecological and evolutionary impacts that 
extend well beyond just the loss of taxonomic diversity (7-9). For 
example, megaherbivore extinctions in the Late Quaternary altered 
plant-herbivore interactions in a way that led to marked shifts in 
plant communities, which, in turn, altered fire regimes and biogeo- 
chemical cycling (9, 10). Herbivores also produce a substantial amount 
of the greenhouse gas methane, and studies have suggested that the 
resulting methane reduction from megaherbivore extinctions in the 
Late Quaternary could have been responsible for a 0.08° to 0.20°C 
decrease in Earth’s climate, leading into the Younger Dryas (8). Al- 
though there is strong evidence that current species loss continues 
to affect large vertebrates disproportionately (11, 12), our under- 
standing of which trophic group is at the greatest risk of extinction 
is primarily driven by anecdotal evidence and correlations with spe- 
cies traits linked to extinction risk. 
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Trophic group, an organism’s position in a food chain, is hypoth- 
esized to be associated with modern-day extinction risk (13-15). Spe- 
cifically, higher trophic groups (e.g., predators) are thought to be 
more at risk than lower ones (e.g., herbivores) (13-16). Extinction 
risk is predicted to correlate with trophic position for several rea- 
sons, including the greater energetic demands (14) and dependence 
of predators on lower trophic levels for food, which themselves may 
be in decline (15, 17). Carbone and Gittleman (18) estimated that 
10,000 kg of prey is needed to support just 90 kg of a given carnivore 
species, and prey depletion has been linked to declines of many large 
carnivores, including tigers, dholes, and several species of leopard 
(17). Declining populations of natural prey combined with retracting 
home ranges also increase human-predator conflicts as predators 
expand into human-occupied areas (19). The tenuous relationship 
between humans and predators is exemplified by the fact that hunt- 
ing and trapping has, at least in part, been implicated in the decline 
of 80% of threatened species in the order Carnivora (1). 

Understanding how risk varies among different trophic groups 
is important because, regardless of the traits driving extinction risk, 
nonrandom patterns of extinction across trophic groups can result 
in changes to trophic structure that may influence many aspects of an 
ecosystem (13). For example, declines in scavenging vultures in India 
have led to increased risk of the spread of diseases such as rabies and 
anthrax (13). In another example, declines in large-bodied frugivores 
in tropical forests are reshaping tree communities, which, in turn, is 
predicted to lead to an overall reduction in forest carbon storage (20). 

Extinction risk is shaped by historical and contemporary socio- 
ecological factors that have affected a region (e.g., time since human 
colonization, gross national product, and industrial development) 
(2, 16). The spatial variability in the human footprint (21) may ulti- 
mately affect trophic groups differently, leading to geographic- or 
habitat-level differences in extinction risk among predators, omni- 
vores, and herbivores. Furthermore, humans have played a role in 
species extinctions since at least the late Pleistocene (4), and patterns 
of threatened species across trophic groups may be derived from past 
extinctions. These processes can cause hysteresis in extinction threat: 
If humans heavily targeted predators in the past, then the predatory 
species that survived are likely to be the most resistant species to 
anthropogenic impacts and thus may be less likely to be threatened 
by contemporary drivers. 
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We examined possible differences in extinction risk across tro- 
phic groups using three different approaches. First, we examined 
contemporary patterns in extinction risk among herbivores, omni- 
vores, and predators. To do this, we focused on comparing threat 
patterns in trophic groups and diets across major taxonomic groups 
(mammals, birds, and reptiles) globally, across different terrestrial 
and oceanic regions, and across different habitat types. Second, we 
examined the role that past extinctions may have played in the trajec- 
tory of contemporary patterns by examining the proportions of re- 
cently extinct mammal, bird, and reptile species and late Pleistocene 
extinct mammals in each trophic group. Last, we determined whether 
trophic group was an important factor driving extinction risk by ex- 
amining how body size and trophic group interact to affect a species’ 
threat status. To examine our objectives, we assembled a database 
categorizing the trophic group (predator, herbivore, and omnivore) 
and diet of all non-data-deficient mammals, birds, and reptiles 
(22,166 species) assessed by the International Union for Conserva- 
tion of Nature (IUCN) (1), the global authority on species extinc- 
tion and extinction risk. 


RESULTS 

Patterns of extinction risk across trophic groups and diets 
The background fraction of all IUCN threatened mammal, bird, 
and reptile species that we examined, calculated by dividing the to- 
tal number of threatened species by the total number of all species, 
was estimated to be 18%. Within these species, we found that herbi- 
vores were the most at-risk trophic group globally, with ~25% [95% 
confidence interval (CI), 24 to 27%] ofall herbivore species listed as 
threatened (Fig. 1A). This amounts to ~300 more herbivores listed 
as threatened than would be expected given the background threat- 
ened fraction of 18%. The proportion of herbivores at risk of extinc- 
tion was higher than omnivores and predators (table S1), with 8% 
(+0.57%) more herbivores listed as threatened compared to omni- 
vores and 10% (+0.52%) more than predators. The pattern of higher 
extinction risk in herbivores relative to other trophic groups was 
consistent even when the three taxonomic classes were analyzed sepa- 
rately (table S1 and Fig. 1A). 

Overall, herbivorous reptiles (e.g., tortoises) had the highest pro- 
portion of at-risk species, with 52% (45 to 59%) listed as threatened. 
It should be noted that as of 2019, ~70% of described reptile species 
have been evaluated by the IUCN, whereas complete assessments 
have been made for mammals and birds. If this subsampling of rep- 
tiles was, in any way, trophically biased, then it could have intro- 
duced bias into our results for reptiles. 

To investigate the potential ecosystem-level effects of species loss 
across trophic groups, we reran our analyses using more detailed 
diet groups. In terms of herbivorous diets, we found that consumers 
of fruit and general plant parts (e.g., leaves, roots, and stems) all had 
higher proportions of threatened species compared to the background 
fraction of 18% (Fig. 1B). Predatory diets that exhibited elevated risk 
of extinction included piscivores and scavengers (Fig. 1B). 


Patterns of extinction risk in trophic groups across 
geographic regions 

To identify region-specific patterns in the extinction risk of different 
trophic groups, we examined the proportions of threatened species 
in herbivores, omnivores, and predators for the five major marine 
regions and 13 land regions classified by the IUCN. The most con- 
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Fig. 1. Patterns of extinction risk by trophic group and diet. (A) Mean propor- 
tions of threatened vertebrate species (+95% CI) within herbivores, omnivores, and 
predators summarized for all species combined and within mammals, birds, and 
reptiles. Letters indicate the results for comparisons across trophic group using 
Tukey's post hoc tests. Differing letters indicate significant differences among trophic 
groups (P < 0.05). (B) Mean proportions of threatened species (95% Cl) within 
each diet group across all species combined. Some species can be classified in mul- 
tiple diet groups. Asterisks indicate where the proportions of threatened species 
are significantly different from the background fraction (dashed horizontal line). 
Numbers indicate sample size. 


sistent pattern shared across geographic regions for all species com- 
bined and for individual taxonomic classes was again an elevated 
vulnerability of herbivores to extinction risk (table S2 and Fig. 2). 
We identified elevated risk (proportion higher than the background 
fraction) for herbivores in at least one taxonomic class in 80% of 
marine regions and 85% of land regions (table S2 and Fig. 2). Her- 
bivorous reptiles, followed by herbivorous mammals, contributed 
the majority of this risk. The second most common at-risk trophic 
group-taxonomic class pairing occurred in omnivorous reptiles, 
which had elevated risk in 80% of marine regions and 46% of land 
regions (table S2 and Fig. 2). Last, predators had elevated levels of 
risk in only 23% of land regions and 40% of the marine regions. Pred- 
atory reptiles and predatory birds (e.g., seabirds) drove the geographic 
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Symbols legend 
a Herbivorous mammals @& Herbivorous reptiles KaD) Omnivorous reptiles 
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YE Predatory birds 
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Fig. 2. Trophic group patterns in extinction risk across land and marine regions. Silhouettes indicate trophic groups and taxonomic classes with higher proportions 
of threatened species in a given region compared to the background fraction. Regions without symbols had either similar or lower proportions of threatened species 
compared to the background fraction. Colors illustrate land and marine regional boundaries. Summary statistics can be found in table S1. 


patterns observed in at-risk predators. Overall, sub-Saharan Africa 
and the Atlantic Ocean had the highest numbers of taxonomic classes 
and trophic groups classified as threatened (Fig. 2). 


Patterns of extinction risk in trophic groups across habitats 
We also examined trophic group patterns in extinction risk across 
five aquatic (marine and inland wetland habitats) and seven terres- 
trial habitats. Consistent with our previous findings, herbivores, fol- 
lowed by omnivores, had the greatest representation of at-risk species 
across habitat types. Herbivores had elevated risk in 100% of aquatic 
habitats and 57% of terrestrial habitats (table S3 and Fig. 3). Herbiv- 
orous reptiles inhabiting aquatic ecosystems were highly threatened, 
with 100% of herbivorous reptiles in marine oceanic, marine intertidal, 
and marine neritic habitats listed as threatened (Fig. 3). In terrestrial 
systems, forests had the highest risk of extinction in herbivores, with 
all three taxonomic classes expressing elevated risk of extinction 
(Fig. 3). Similar to herbivores, omnivores also showed high levels of 
risk in all aquatic habitats. Patterns of extinction risk across habitats 
for omnivores were driven almost entirely by reptiles, except in ma- 
rine neritic habitats where both omnivorous reptiles and mammals 
showed elevated risks of extinction (table S3 and Fig. 3). Higher pro- 
portions of threatened predators occurred in 60% of aquatic habi- 
tats and 14% of terrestrial habitats (table S3 and Fig. 4). Predatory 
birds, especially those in marine habitats (e.g., seabirds), drove our 
observed patterns in extinction risk in predators (table S3 and Fig. 3). 


Historical patterns in trophic group extinctions 

To examine whether and how past extinctions have influenced cur- 
rent patterns in threatened species, we compared the proportions of 
recently extinct mammal, bird, and reptile species (i-e., species clas- 
sified as extinct or extinct in the wild in the last 500 years by the 
IUCN) and late Pleistocene extinct mammals (beginning 11,000 years 
ago for Africa, North America, and South America and 50,000 years 
ago for Australia) in each of our three trophic groups (22). In both 
recently extinct species and late Pleistocene extinct mammals, we 
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found herbivores to have the greatest proportion of extinct species 
(table S4 and Fig. 4). In contrast, predators had lower proportions 
of extinctions compared to herbivores and, except for recently ex- 
tinct birds, had lower levels of extinction compared to the background 
proportion. Omnivores had similar levels of extinction compared to 
the background proportion, with the exception of Pleistocene ex- 
tinct mammals, which had a significantly lower proportion (table S4 
and Fig. 4). 


Potential drivers of trophic group patterns in extinction risk 
Our above analyses show trophic group patterns in extinction risk 
across multiple spatial and temporal scales, but it does not deter- 
mine whether trophic group itself is an important factor driving ex- 
tinction risk. To answer this question, we ran phylogenetic generalized 
linear models (GLMs) to test associations between extinction risk, 
body size, and trophic group. We found evidence that trophic group 
and body size are important factors driving the extinction risk pat- 
terns observed in our study. However, the overall importance and 
relationship between these two factors varied for mammals, birds, 
and reptiles (Fig. 5 and table S5). 

For mammals, we found that body mass and trophic group were 
both important factors that independently affected threat status, with 
large-bodied organisms and herbivores experiencing the greatest risk 
of extinction (table S5 and Fig. 5). Neither individual anthropogenic 
threat drivers (e.g., biological resource use and habitat alteration) 
nor the total number of anthropogenic drivers experienced by a tro- 
phic group helped explain the overrepresentation of threatened her- 
bivores in mammals (table S6 and Fig. 6). In many cases, herbivorous 
mammals were the least affected or similarly affected by a driver com- 
pared to predatory and omnivorous mammals (table S6 and Fig. 6). 

In birds, we found that body mass was positively correlated with 
threat status. However, trophic group only marginally affected the 
threat status of birds, with herbivores having a slightly higher pro- 
portion of threatened species compared to predators (table S5 and 
Fig. 5). In birds, small-bodied herbivores were more affected by habitat 
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Fig. 3. Trophic group patterns in extinction risk across habitat types. Icons indicate trophic groups with lower proportions (white), similar proportions (gray), or 
higher proportions (black) of threatened species compared to the background fraction. Numbers indicate the mean proportion (+95% Cl) listed as threatened for each 
trophic group-habitat combination. Letters indicate the results for comparisons across trophic group using Tukey's post hoc tests. Differing letters indicate significant 
differences among trophic groups (P< 0.05). Summary statistics can be found in table S2. 


alteration and pollution compared to omnivores (Fig. 6). Small her- 
bivorous birds were also more affected by invasive species and pol- 
lution compared to small predatory birds; however, this pattern was 
the reverse in large-bodied birds (table S7 and Fig. 6). 

In reptiles, trophic group mediated the association between body 
mass and threat status. Here, we found that the probability of being 
threatened was positively related to body mass in herbivores and 
omnivores, but that the inverse was true for predatory reptiles (table 
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S5 and Fig. 5). For reptiles, invasive species was the only anthropogenic 
pressure to significantly affect herbivores more than other trophic 
groups (table S8 and Fig. 6). 


DISCUSSION 
Numerous studies have helped to identify species’ traits that cor- 
relate with extinction risk to shed light on current patterns in the 
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Fig. 4. Historical patterns in extinctions across trophic groups. The mean pro- 
portions (+95% CI) of extinct species for recently extinct birds, mammals, and reptiles 
combined (all species) and separately, and Pleistocene extinct mammals. Asterisks 
indicate where the proportions of extinct species are significantly different from 
the background fraction (dashed horizontal line). Numbers indicate total number 
of species within each group. Letters indicate the results for comparisons across 
trophic group using Tukey's post hoc tests. Differing letters indicate significant dif- 
ferences among trophic groups (P< 0.05). 


Anthropocene extinction crisis (12, 14, 15). Findings from some of 
these past studies, alongside a host of single species research fo- 
cused on extinction threat in charismatic predators, have led to the 
general assumption that predators are, for a variety of reasons, at a 


Atwood etal., Sci. Adv. 2020; 6: eabb8458 5 August 2020 


oO 
=] — Herbivores 
— Omnivores 
2 4 — Predators 
o w 
ei) 3 
~ | 3 
2 D. 
v v 
a4 
oO 
8 AF 
0 2 4 6 8 
To 
O C4 
co 
Lox 
O°] 
L o w 
59| = 
cx z 
O°] 
= 
“Ln 
Oo] 
Q 
Oc] 
a T T T T T 
0 1 2 3 4 5 
aj 
co 
a7 
© | J 
S D 
ze} 
ae a 
oO — 
oO 
S Wn 
a4 
oO 
o Hr 


-2 0 2 4 6 
Fig. 5. Relationships between body mass, trophic group, and threat status. 
The relationship between body size, trophic group, and the proportion of mam- 
mals, birds, and reptiles listed as threatened by the IUCN. Models account for phy- 
logenetic relatedness. Solid lines represent means, while shaded areas represent 
95% Cl. See table S5 for summary statistics. 


higher risk of extinction than other trophic groups (14-16, 19). By 
systematically examining the patterns of at-risk species across dif- 
ferent trophic groups, our study identified that the threat of extinction 
is trophically skewed. However, contrary to many of these previous 
expectations, we found herbivores to be the most at-risk trophic group 
among mammals, birds, and reptiles. Although geography and habitat 
influenced taxonomic class-specific results, herbivores consistently 
had the highest representation of at-risk species in the present day, the 
recent past, and the late Pleistocene. In many cases, reptiles were the 
primary group driving our observed patterns in extinction risk of 
modern herbivores. This result is noteworthy, given that much has yet 
to be determined regarding their contemporary functional roles in 
ecosystems (23). 

We identified a few instances where predators show elevated risk 
of extinction. When more detailed diet categories were analyzed, 
both piscivores and scavengers had elevated risk compared to back- 
ground levels. Consistent with the higher extinction risk in piscivores, 
we found that elevated risk of extinction in predators occurred al- 
most exclusively in marine habitats, which suggests that extinction 
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Fig. 6. Impacts of anthropogenic drivers of global change on threatened species. The relationship between body size, trophic group, and the proportion of species 
affected by habitat alteration, biological resource use, climate change, pollution, and invasive species. Models account for phylogenetic relatedness. Solid lines represent 
means, while shaded areas represent 95% CI. See tables S6 to S8 for summary statistics. 


pressures threatening predators may be greater in the ocean than on 
land. Past research has proposed that, in many regions, humans are 
preferentially and unsustainably exploiting marine organisms at the 
top of the food chain (24). However, this research has focused mainly 
on fish (24), which we did not include in our study. A likely produc- 
tive future research endeavor would be to determine how the inclu- 
sion of bony and cartilaginous fishes and other taxonomic classes not 
in this study (e.g., invertebrates) would influence interpretations about 
how trophic group shapes extinction risk overall, as well as how it 
influences our observations about the relative risk for different tro- 
phic groups in terrestrial and aquatic environments. 

Although we were unable to identify a single anthropogenic driver 
for the global decline in herbivores, we did find that certain drivers 
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disproportionately affect some groups of herbivores more relative 
to other trophic groups. For example, we found that invasive spe- 
cies affect herbivorous reptiles disproportionately compared to om- 
nivores and predators. Invasive vertebrates (e.g., rats), insects (e.g., 
fire ants), and plants (e.g., Hottentot fig) have all been implicated in 
the decline and even extinction of several reptiles (25). Furthermore, 
we found that invasive species, pollution, and habitat alteration af- 
fect small herbivorous birds disproportionately. However, for both 
of these cases, it is unclear why these anthropogenic drivers would 
target herbivores more than other trophic groups. A new challenge 
for conservation biology will be to identify clear mechanisms respon- 
sible for this apparent association between trophic group and ex- 
tinction risk. These studies are likely to include investigations into 
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the interactive effects of multiple anthropogenic drivers and their im- 
pact upon intrinsic species traits that are associated with herbivory. 

Our study determined that trophic group is an important factor 
driving extinction risk in reptiles, mammals, and birds. Past studies 
examining the effect of trophic group or diet on extinction risk have 
found either no effect of trophic group (26) or that higher-order pred- 
ators are at the highest risk of extinction (14, 16). These past studies, 
however, focused on select groups of organisms like animals in the 
order Carnivora (16), or squamate reptiles (26), suggesting that these 
more focused studies mask the more general role of trophic group 
as a significant variable that influences extinction patterns in mam- 
mals, birds, and reptiles collectively. Critically, our results show that 
any future endeavors to explain patterns of extinction risk in mam- 
mals, birds, and reptiles need to account for trophic group. 

In addition to trophic group, our study found that body size is an 
important trait that independently drives extinction risk in mam- 
mals and birds and interacts with trophic group to influence extinc- 
tion risk in reptiles. Since the late Pleistocene, the average body mass 
of organisms has been declining because of size-selective threats that 
are robustly linked to human activities (27). Our results agree with 
findings from past studies that the selection against large-bodied or- 
ganisms is likely to continue under business as usual management 
scenarios (11, 12). Our study adds to this body of literature by show- 
ing that, except for predatory reptiles, large-bodied organisms across 
all three trophic groups are at a high risk of extinction, but herbi- 
vores are disproportionately the most at-risk trophic group within 
these large-bodied species. The ecological effects of the loss of large- 
bodied herbivores will depend, to some degree, on the capacity of 
smaller-bodied species to compensate numerically and/or functionally 
(28). Many body size-specific functions (e.g., dispersal of large-seeded 
fruits by sloths) cannot, however, be replicated by simply increasing 
abundances of smaller-bodied species (28). 

Prehistoric extinctions of megaherbivores drastically changed the 
structure and functioning of Earth’s ecosystems by altering vegeta- 
tion dynamics, fire regimes, carbon cycling, and biogeochemical cy- 
cling (7-9). Our results highlight that Earth is once again experiencing 
declines and extinctions that are disproportionately affecting large 
herbivores. However, how these declines and extinctions of herbi- 
vores are likely to affect the trajectory of life on Earth is not yet known, 
but studies have linked modern herbivores with ecosystem process- 
es as diverse as the evolution of plant and predator traits, ecosystem 
resilience/resistance, nutrient cycling, fire regimes, greenhouse gas 
dynamics, plant regeneration, and primary production (29). 


MATERIALS AND METHODS 

We created a database containing biological, ecological, and geo- 
graphical data on all extant mammal (4858), bird (10,910), and reptile 
(6398) species assessed on the 2019-2 IUCN Red List of Threatened 
Species (1) as well as recently extinct (extinct in the wild and extinct on the 
IUCN Red List of Threatened Species) mammals, birds, and reptiles 
and late Pleistocene extinct mammals derived from McCauley et al. 
(2) and Smith et al. (27). First, we assigned binary threat assessments 
to each species. Second, we characterized the diets of each species 
based on available information. Third, we classified the trophic groups 
of each species based on our diet characterization. We then com- 
bined these data with species-specific geographic location, habitat, 
and threat information from the IUCN Red List of Threatened Spe- 
cies. We used the completed database to quantify the observed pro- 
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portions of threatened or extinct species in specific trophic groups 
(i.e., predator, omnivore, and herbivore) in relation to each other 
and the background fractions of threatened or extinct species. 


Binary threat classification 

For each extant species assessed by the IUCN, we assigned a binary 
threat classification of “threatened” or “nonthreatened.” Species with 
IUCN Red List assessments of “critically endangered,” “endangered,” 
and “vulnerable” were considered threatened with extinction, and 
species with “near threatened,” “least concern,” and “lower risk” as- 
sessments were considered nonthreatened. This binary classification 
system represents the two potential outcomes in an extinction event 
(i.e., extinction or survival) and is consistent with analytical approaches 
previously used for examining extinction risk (11). Species assessed 
as “extinct” or “extinct in the wild” were not included in threatened 
status analyses. Instead, they were classified as recently extinct in 
analyses looking at historical trophic group patterns. Species cate- 
gorized as data deficient by IUCN were not used in our analyses. 
Mammals from the extinct species lists from McCauley et al. (2) and 
Smith et al. (27) were classified as extinct in the late Pleistocene and 
were only used in analyses of historical trophic group patterns. 


Diet categories 

We classified the diet of all mammal, bird, and reptile species assessed 
by the IUCN using published literature, reference texts, databases, or 
extrapolation from related taxonomic class. Species’ diets were de- 
termined from various data types (i.e., percent, ranked presence- 
absence, and unranked presence-absence of specific food items; table 
S9). We classified diet categories using a binary response of present 
or absent. Only adult diet information from wild populations was 
used. The primary diet categories (with subcategories in parentheses) 
were endotherm (mammal and bird), herptile (amphibian and rep- 
tile), fish, invertebrate (insect), carrion, plant (foliage, root, and wood), 
fruit, nectar, and grain (30, 31). Carrion eaters (i.e., scavengers) were 
defined as species that consume animal-based carrion. To be classi- 
fied in the plant category, a species must consume algae, fungus, leaves, 
shoots, roots, wood, flowers, pollen, sepals, or other miscellaneous 
plant material (e.g., vegetable matter) (30). 

For diets described in the literature as percentages, we used spe- 
cific percentage cutoffs to classify diet categories as either present or 
absent. Percentage cutoffs were adjusted on the basis of the number 
of unique diet items included in a species’ diet, because as the num- 
ber of diet items increases, the weight of a diet category increases. 
For example, if a species only consumed two diet items, and one 
diet item comprised <20% of the diet, then that item was classified 
as absent. However, if the species’ diet had more than two diet items, 
then any item with a >20% contribution was classified as present. 

For diets described as presence-absence data, diet categories were 
classified according to keywords. We included diet items associated 
with the following keywords: “primary,” “secondary,” “mostly,” “also,” 
“frequently,” “regularly,” “usually,” and “fair amount.” We excluded 
diet items associated with the keywords “tertiary,” “sometimes,” “oc- 
casionally,” “rarely,” “small quantities,” “and even,” “opportunisti- 
cally,” and “at times supplemented with.” If a diet item was associated 
with uncertainty statements of “probably,” “possibly,” “may,” or “pre- 
sumably,” then it was excluded unless no other diet information was 
available. For diets described using unranked presence-absence data 
(e.g., “eats fruit, insects, and seeds”), all listed diet items were classi- 
fied as present. 
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Trophic group 

Using dietary information, we classified species into three trophic 
groups: predator, omnivore, and herbivore. For ranked presence- 
absence data, predators were defined as animals that primarily con- 
sumed animal-based diet items and could consume plant-based diet 
items in small amounts; herbivores primarily consumed plant-based 
diet items and could consume animal-based diet items in small 
amounts, and omnivores primarily or secondarily consumed a mix- 
ture of both types of diet items. For unranked presence-absence data, 
predators consumed only animal-based diet items, herbivores con- 
sumed only plant-based diet items, and omnivore diets included 
items from both categories. To explore the sensitivity of results to 
our classification of trophic groups, trophic groups were classified 
in three different ways where the contributions of the herbivore or 
predator diet shifted: 70% cutoff, 80% cutoff, and 90% cutoff. Results 
described in the main text represent the most conservative view of 
predators and herbivores, with predators defined as species that con- 
sume a >90% animal-based diet, herbivores as species that consume 
a >90% plant-based diet, and omnivores as species that consume an 
11 to 89% animal or plant-based diet. The 80% diet defines preda- 
tors as species that consume a >80% animal-based diet, herbivores 
as species that consume a >80% plant-based diet, and omnivores as 
species that consume a 21 to 79% animal or plant-based diet. The 
70% diet defines predators as species that consume a >70% animal- 
based diet, herbivores as species that consume a >70% plant-based 
diet, and omnivores as species that consume a 31 to 69% animal or 
plant-based diet. In general, our conclusions about how trophic group 
shaped extinction risk were insensitive to our methods for defining 
trophic groups (fig. S1). 

When diet information could not be found for a species, we ex- 
trapolated trophic group from the sister taxonomic class. If a spe- 
cies was formerly classified as another species or a subspecies, then 
the trophic group from the former classification was used. If the species 
was not formally classified as another species or a subspecies, then 
trophic group was extrapolated from the closest living congener 
or confamiliar using previously published phylogenies (32-34). To 
ground-truth this approach, we first removed all extrapolated spe- 
cies from our dataset and then split the remaining data into training 
data (75% of the data) and test data (25%). We then assigned tro- 
phic group to species in our test data set using the trophic group of 
the closest living relative in our training data. We repeated this pro- 
cess across 1000 bootstrapped phylogenies to provide a mean accu- 
racy (+SE) and compared our extrapolated trophic group results 
with the known trophic groups in the original data set. We found 
that our accuracies for predicting the trophic group of mammals, 
birds, and reptiles were 94% (+0.02), 82% (+0.03), and 88% (+0.05), 
respectively. 

For late Pleistocene extinct mammals as listed by Smith et al. (27), 
only trophic group was classified and not diet due to a paucity of 
specific diet information. Trophic group was categorized on the basis 
information collated from the primary literature, databases, and closest 
extant congeners using similar protocols to the dietary categorization 
described previously. 


Regions, habitats, and anthropogenic global change drives 

To explore anthropogenic effects on trophic group patterns and the 
variability in threat classifications, we downloaded all species-specific 
extrinsic variables (i.e., habitats, regions, and associated anthropo- 
genic threats) from the IUCN Red List of Threatened Species. Major 
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geographic (land and marine) regions, habitats, and anthropogenic 
threats were classified as either present or absent for each species. 
The major anthropogenic global change drives classified by the IUCN 
were simplified into five major anthropogenic threat categories (3): 
habitat alteration (residential and commercial development, agri- 
culture and aquaculture, energy production and mining, transpor- 
tation and service corridors, human intrusions and disturbance, or 
natural system modifications), biological resource use (e.g., overex- 
ploitation), climate change, pollution, and invasive species. 


Body size 

We used adult body size for classification. When male and female 
body sizes were given, we used the geometric mean of the two val- 
ues. When a range of body sizes was given, we used the geometric 
mean of the smallest and largest value (30). Body mass was most 
readily available for birds and mammals. However, body size mea- 
surements for reptiles were not standardized across taxonomic 
orders. Body size for Squamata was commonly reported as total length 
or snout-vent length, while body size for Testudines was reported as 
standard carapace length. To overcome nonstandardization methods 
for measuring body size, we converted all length measurements to 
body mass using previously published family-level length-mass re- 
gressions. References for length-mass regressions are provided in the 
Supplementary Materials. If previously published regressions were 
not available, then we developed regressions from species in our 
data that contained both body mass and length data. For species 
with missing body size data, we extrapolated body size using the 
average of all congeners (excluding congeners that had been previ- 
ously extrapolated this way). 


Analyses 

The IUCN Red List of Threatened Species consists of imperfect data 
as many species are classified as data deficient, and others have not 
been described at all. Because of the imperfect nature of the IUCN 
data, we used the available data to generate mean estimates and 95% 
confidence limits of the proportions of threatened species for differ- 
ent trophic groups. Therefore, although we may not know what the 
exact true proportion of threatened species is for each trophic group, 
our CIs indicate boundaries within which the true number likely 
lies. We investigated patterns in the proportion of threatened spe- 
cies for each trophic group (herbivore, omnivore, and predator) at 
the global scale for extant vertebrate species (mammals, birds, and 
reptiles), recently extinct vertebrate species (mammals, birds, and 
reptiles), and Pleistocene extinct vertebrate species (mammals only). 
Differences in proportions of threatened species among trophic groups 
were analyzed using GLMs, with threat status as the response vari- 
able and trophic group as the predictor variable. We used a binomial 
error structure (threatened or not threatened) and a logit link func- 
tion. Each species was treated as a single replicate and recorded as 
either threatened or not, or when historical patterns were examined, 
extinct or not. Following the GLM, post hoc Tukey’s tests were used 
to identify differences among trophic groups in the proportions of 
threatened or extinct species. For each trophic group, we produced 
mean estimates and 95% CIs of the proportion of threatened species 
using the “predict” function in R. These 95% CIs were compared to 
the proportion of threatened or extinct species across background 
levels, which were calculated by taking the number of threatened spe- 
cies and dividing it by the total number of species in that analysis. If 
the 95% CI did not intersect the background fraction of threatened 
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species (i.e., background level), then the trophic group was desig- 
nated as having a higher or lower proportion of threatened species 
than the background fraction. After running the analysis across all 
taxa in the IUCN Red List of Threatened Species, we split the data 
on the basis of taxonomic class and ran individual models for each 
class separately. 

To understand the relationships between specific diet types and 
the proportion of threatened species, we ran subsequent analyses 
where each diet type was treated as a binary descriptive variable. 
This technique meant that each species could be classified in multi- 
ple diet categories (i.e., an herbivore could be classed as consuming 
nectar and grain), meaning that the diet categories were not mutually 
exclusive. We used a GLM to understand how different diet catego- 
ries are compared to the background fraction of threatened species, 
and again used the predict function to generate means and 95% CIs 
for each diet type. These 95% CIs for each specific diet were then 
compared to the background fraction of threatened species. 

Following the global analyses, we split our data on the basis of 
habitat type and geographic region and repeated the analyses to un- 
derstand how patterns in threat status of different trophic groups 
varied across geographic regions and habitat types. For each habitat 
or geographic region, we repeated the techniques described above 
for the global analyses. First, we used a GLM and post hoc Tukey’s 
test to assess differences in the proportion of threatened species among 
trophic groups for each habitat or region. We then used the predict 
function to generate mean proportions and 95% CIs for each taxo- 
nomic class and compared them to the proportion of threatened spe- 
cies within that taxonomic class (either all species, mammals, birds, 
or reptiles) across the globe. All analyses were carried out in R. 

To determine potential drivers of extinction risk, we ran phylo- 
genetic GLMs with both trophic group and body size as predictor 
variables. To determine whether specific anthropogenic drivers or 
the total number of anthropogenic drivers disproportionately affected 
herbivores, we examined the proportion of threatened species in each 
trophic group affected by resource use, climate change, habitat al- 
teration, invasive species, and pollution while accounting for both 
body mass and phylogeny. Phylogenetic signal is well documented 
in IUCN threat status data (33, 35), and to ensure that this would 
not affect our models assessing support for potential mechanisms (36), 
we used a phylogenetic comparative approach. Our use of phyloge- 
netic GLMs (37), as implemented in phylolm (38), ensures that our 
model estimates are robust to species shared evolutionary history. 
Further, to account for uncertainty about phylogenetic topology or 
divergence dates, we repeated our analyses across Bayesian posterior 
distributions of phylogenetic trees [reviewed in (39)]. We took 1000 
phylogenies for each taxonomic group—Faurby and Svenning (32) 
for mammals, Jetz et al. (34) for birds, and Tonini et al. (33) for 
reptiles—and repeated all analyses across them, reporting mean pa- 
rameter estimates across all these models. To ensure that our trait 
and phylogenetic data were compatible, we matched all phylogenetic 
data to the taxonomy used by IUCN (that is, the basis of our trait 
data; see above) using taxize (40). We used Akaike information cri- 
terion to select the most parsimonious model. R codes to repeat our 
analyses are given in the Supplementary Materials. 


SUPPLEMENTARY MATERIALS 


Supplementary material for this article is available at http://advances.sciencemag.org/cgi/ 
content/full/6/32/eabb8458/DC1 


Atwood etal., Sci. Adv. 2020; 6: eabb8458 5 August 2020 


REFERENCES AND NOTES 

1. IUCN, The IUCN Red List of Threatened Species. Version 2019-2 (IUCN, 2019); 
http://www.iucnredlist.org. 

2. D.J. McCauley, M. L. Pinsky, S. R. Palumbi, J. A. Estes, F. H. Joyce, R. R. Warner, Marine 
defaunation: Animal loss in the global ocean. Science 347, 1255641 (2015). 

3. H.S. Young, D. J. Mccauley, M. Galetti, R. Dirzo, Patterns, causes, and consequences 
of Anthropocene defaunation. Annu. Rev. Ecol. Evol. Syst. 47, 333-358 (2016). 

4. A.D.Barnosky, P. L. Koch, R. S. Feranec, S. L. Wing, A. B. Shabel, Assessing the causes 
of late Pleistocene extinctions on the continents. Science 306, 70-75 (2004). 

5. M. Urban, Accelerating extinction risk from climate change. Science 348, 571-573 
(2015). 

6. D.U. Hooper, E. C. Adair, B. J. Cardinale, J. E. K. Byrnes, B. A. Hungate, K. L. Matulich, 

A. Gonzalez, J. E. Duffy, L. Gamfeldt, M. I. O'Connor, A global synthesis reveals biodiversity 
loss as a major driver of ecosystem change. Nature 486, 105-108 (2012). 

7. A. D. Barnosky, E. L. Lindsey, N. A. Villavicencio, E. Bostelmann, E. A. Hadly, J. Wanket, 

C. R. Marshall, Variable impact of late-Quaternary megafaunal extinction in causing 
ecological state shifts in North and South America. Proc. Natl. Acad. Sci. U.S.A. 113, 
856-861 (2016). 

8. F.A. Smith, J. I. Hammond, M. A. Balk, S. M. Elliott, S. K. Lyons, M. I. Pardi, C. P. Tomé, 

P. J. Wagner, M. L. Westover, Exploring the influence of ancient and historic 
megaherbivore extirpations on the global methane budget. Proc. Natl. Acad. Sci. U.S.A. 
113, 874-879 (2016). 

9. J.L. Gill, Ecological impacts of the late Quaternary megaherbivore extinctions. 

New Phytol. 201, 1163-1169 (2014). 

0. J.L. Gill, J. W. Williams, S. T. Jackson, K. B. Lininger, G. S. Robinson, Pleistocene 
megafaunal collapse, novel plant communities, and enhanced fire regimes in North 
America. Science 326, 1100-1104 (2009). 

1. W. J. Ripple, C. Wolf, T. M. Newsome, M. Hoffmann, A. J. Wirsing, D. J. McCauley, Extinction 
risk is most acute for the world’s largest and smallest vertebrates. Proc. Natl. Acad. Sci. 
U.S.A. 114, 10678-10683 (2017). 

2. M. Cardillo, G. M. Mace, K. E. Jones, J. Bielby, O. R. P. Bininda-Emonds, W. Sechrest, 

C. D. L. Orme, A. Purvis, Multiple causes of high extinction risk in large mammal species. 
Science 309, 1239-1241 (2005). 
3. J.A. Estes, J. Terborgh, J. S. Brashares, M. E. Power, J. Berger, W. J. Bond, S. R. Carpenter, 
T. E. Essington, R. D. Holt, J. B. C. Jackson, R. J. Marquis, L. Oksanen, T. Oksanen, R. T. Paine, 
E. K. Pikitch, W. J. Ripple, S. A. Sandin, M. Scheffer, T. W. Schoener, J. B. Shurin, 
A. R. E. Sinclair, M. E. Soulé, R. Virtanen, D. A. Wardle, Trophic downgrading of planet 
Earth. Science 333, 301-306 (2011). 
14. A. Purvis, J. L. Gittleman, G. Cowlishaw, G. M. Mace, Predicting extinction risk in declining 
species. Proc. Biol. Sci. 267, 1947-1952 (2000). 

15. R. D. Holt, J. H. Lawton, G. A. Polis, N. D. Martinez, Trophic rank and the species-area 
relationship. Ecology 80, 1495-1504 (1999). 

16. M. Cardillo, A. Purvis, W. Sechrest, J. L. Gittleman, J. Bielby, G. M. Mace, Human 
population density and extinction risk in the world’s carnivores. PLOS Biol. 2, e197 
(2004). 

17. C. Wolf, W. J. Ripple, Prey depletion as a threat to the world’s large carnivores. R. Soc. Open 
Sci. 3, 160252 (2016). 

18. C. Carbone, J. L. Gittleman, A common rule for the scaling of carnivore density. Science 
295, 2273-2276 (2002). 

19. W.J.Ripple, J. A. Estes, R. L. Beschta, C. C. Wilmers, E. G. Ritchie, M. Hebblewhite, J. Berger, 
B. Elmhagen, M. Letnic, M. P. Nelson, O. J. Schmitz, D. W. Smith, A. D. Wallach, 

A. J. Wirsing, Status and ecological effects of the world’s largest carnivores. Science 343, 
1241484 (2014). 

20. C. Bello, M. Galetti, M. A. Pizo, L. F. S. Magnago, M. F. Rocha, R. A. F. Lima, C. A. Peres, 

O. Ovaskainen, P. Jordano, Defaunation affects carbon storage in tropical forests. Sci. Adv. 
1, e1501105 (2015). 

21. O. Venter, E. W. Sanderson, A. Magrach, J. R. Allan, J. Beher, K. R. Jones, H. P. Possingham, 
W. F. Laurance, P. Wood, B. M. Fekete, M. A. Levy, J. E. M. Watson, Sixteen years of change 
in the global terrestrial human footprint and implications for biodiversity conservation. 
Nat. Commun. 7, 12558 (2016). 

22. F.A. Smith, S. K. Lyons, S. K. M. Ernest, K. E. Jones, D. M. Kaufman, T. Dayan, 

P. A. Marquet, J. H. Brown, J. P. Haskell, Body mass of late Quaternary mammals. Ecology 
84, 3403 (2003). 

23. E.B. P. de Miranda, The plight of reptiles as ecological actors in the tropics. Front. Ecol. 
Evol. 5, 159 (2017). 

24. R. A. Myers, B. Worm, Rapid worldwide depletion of predatory fish communities. Nature 
423, 280-283 (2003). 

25. J. W. Gibbons, D. E. Scott, T. J. Ryan, K. A. Buhlmann, T. D. Tuberville, B. S. Metts, 

J. L. Greene, T. Mills, Y. Leiden, S. Poppy, C. T. Winne, The global decline of reptiles, 
Déjà Vu Amphibians: Reptile species are declining on a global scale. Six significant threats 
to reptile populations are habitat loss and degradation, introduced invasive species, 


9of10 


ozoz ‘Le IsnBny uo /B10'Bewasuans'sooueape//:dyy wo papeojumogq 


SCIENCE ADVANCES | RESEARCH ARTICLE 


environmental pollution, disease, unsustainable use, and global climate change. 
BioScience 50, 653-666 (2000). 

26. M. Böhm, R. Williams, H. R. Bramhall, K. M. Mcmillan, A. D. Davidson, A. Garcia, L. M. Bland, 
J. Bielby, B. Collen, Correlates of extinction risk in squamate reptiles: The relative 
importance of biology, geography, threat and range size. Glob. Ecol. Biogeogr. 25, 
391-405 (2016). 

27. F.A. Smith, R. E. E. Smith, S. K. Lyons, J. L. Payne, Body size downgrading of mammals 
over the late Quaternary. Science 313, 310-313 (2018). 

28. N. Sekar, C.-L. Lee, R. Sukumar, Functional nonredundancy of elephants in a disturbed 
tropical forest. Conserv. Biol. 31, 1152-1162 (2017). 

29. E.S. Forbes, J. H. Cushman, D. E. Burkepile, T. P. Young, M. Klope, H. S. Young, 
Synthesizing the effects of large, wild herbivore exclusion on ecosystem function. 
Funct. Ecol. 33, 1597-1610 (2019). 

30. H. Wilman, J. Belmaker, S. Jennifer, C. de la Rosa, M. M. Rivadeneira, W. Jetz, EltonTraits 
1.0: Species-level foraging attributes of the world’s birds and mammals. Ecology 95, 2027 
(2014). 

31. W.D. Kissling, L. Dalby, C. Flajgaard, J. Lenoir, B. Sandel, C. Sandom, K. Trajelsgaard, 

J.-C. Svenning, Establishing macroecological trait datasets: Digitalization, extrapolation, 
and validation of diet preferences in terrestrial mammals worldwide. Ecol. Evol. 4, 
2913-2930 (2014). 

32. S. Faurby, J.-C. Svenning, A species-level phylogeny of all extant and late Quaternary 
extinct mammals using a novel heuristic-hierarchical Bayesian approach. Mol. Phylogenet. 
Evol. 84, 14-26 (2015). 

33. J. F.R. Tonini, K. H. Beard, R. B. Ferreira, W. Jetz, R. A. Pyron, Fully-sampled phylogenies of 
squamates reveal evolutionary patterns in threat status. Biol. Conserv. 204, 23-31 (2016). 

34. W. Jetz, G. H. Thomas, J. B. Joy, K. Hartmann, A. O. Mooers, The global diversity of birds 
in space and time. Nature 491, 444-448 (2012). 

35. S. Cannicci, D. Burrows, S. Fratini, T. J. Smith III, J. Offenberg, F. Dahdouh-Guebas, Faunal 
impact on vegetation structure and ecosystem function in mangrove forests: A review. 
Aquat. Bot. 89, 186-200 (2008). 

36. N. Cooper, W. Jetz, R. P. Freckleton, Phylogenetic comparative approaches for studying 
niche conservatism. J. Evol. Biol. 23, 2529-2539 (2010). 

37. A.R. Ives, T. Garland Jr., Phylogenetic logistic regression for binary dependent variables. 
Syst. Biol. 59, 9-26 (2010). 


Atwood etal., Sci. Adv. 2020; 6: eabb8458 5 August 2020 


38. L.S. Tung Ho, C. Ané, A linear-time algorithm for Gaussian and non-Gaussian trait 
evolution models. Syst. Biol. 63, 397-408 (2014). 

39. J.P. Huelsenbeck, F. Ronquist, R. Nielsen, J. P. Bollback, Bayesian inference of phylogeny 
and its impact on evolutionary biology. Science 294, 2310-2314 (2001). 

40. S.A. Chamberlain, E. Sz6cs, Taxize: Taxonomic search and retrieval in R. F1000Res. 2, 191 
(2013). 


Acknowledgments: Symbols in Figs. 1 to 5 were provided by Freepik, and images in Fig. 6 are 
courtesy of the Integration and Application Network, University of Maryland Center for 
Environmental Science (ian.umces.edu/symbols/). Funding: Funding was provided by a USU 
Ecology Center and Graduate Enhancement Award to S.A.V.; an Early Career Research 
Fellowship from the Gulf Research Program of the National Academies of Sciences, 
Engineering, and Medicine to T.B.A.; a gift to UC Santa Barbara from Marc and Lynne Benioff in 
the form of support for the Benioff Ocean Initiative to D.J.M.; Utah State University-Utah 
Agricultural Experiment Station (journal paper number 9282); and the NSF (ABI-1759965 and 
EF-1802605) and the U.S. Forest Service (18-CS-11046000-041) grants to W.D.P. The content is 
solely the responsibility of the authors and does not necessarily represent the official views of 
the Gulf Research Program of the National Academies of Sciences, Engineering, and Medicine. 
Author contributions: T.B.A. and S.A.V. contributed equally to the work. T.B.A., E.H., E.M.P.M., 
DJ.M., and S.A.V. designed the study; S.A.V. and K.H.B. collected data; E.H. and W.D.P. analyzed 
data; and S.A.V. and T.B.A. wrote the manuscript. All authors contributed to editing of the 
manuscript. Competing interests: The authors declare that they have no competing interests. 
Data and materials availability: All data needed to evaluate the conclusions in the paper 
are present in the paper and/or the Supplementary Materials. Additional data related to this 
paper may be requested from the authors. 


Submitted 22 March 2020 
Accepted 25 June 2020 
Published 5 August 2020 
10.1126/sciadv.abb8458 


Citation: T. B. Atwood, S. A. Valentine, E. Hammill, D. J. McCauley, E. M. P. Madin, K. H. Beard, 


W. D. Pearse, Herbivores at the highest risk of extinction among mammals, birds, and reptiles. 
Sci. Adv. 6, eabb8458 (2020). 


10 of 10 


ozoz ‘LE IsnBny uo /Hio bewseousios ’ssoueApe//:dyjy woi papeojumogq 


Science Advances 


Herbivores at the highest risk of extinction among mammals, birds, and reptiles 


Trisha B. Atwood, Shaley A. Valentine, Edd Hammill, Douglas J. McCauley, Elizabeth M. P. Madin, Karen H. Beard and 
William D. Pearse 


Sci Adv 6 (32), eabb8458. 
DOI: 10.1126/sciadv.abb8458 


ARTICLE TOOLS http://advances.sciencemag.org/content/6/32/eabb8458 
SUPPLEMENTARY http://advances.sciencemag.org/content/suppl/2020/08/03/6.32.eabb8458.DC1 
REFERENCES This article cites 39 articles, 13 of which you can access for free 


http://advances.sciencemag.org/content/6/32/eabb8458#BIBL 


PERMISSIONS http:/www.sciencemag.org/help/reprints-and-permissions 


Use of this article is subject to the Terms of Service 


Science Advances (ISSN 2375-2548) is published by the American Association for the Advancement of Science, 1200 New 
York Avenue NW, Washington, DC 20005. The title Science Advances is a registered trademark of AAAS. 


Copyright © 2020 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of 
Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution NonCommercial 
License 4.0 (CC BY-NC). 


ozoz ‘LE IsnBny uo /Hio bewsousios ’ssoueApe//:dyjy wo papeojumogq 


